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Abstract 

The possibility for the occurrence in crystals of a phenomenon, resembling turbu- 
lence, is discussed. This phenomenon, called heterophase turbulence, is manifested by 
the fluctuational appearance inside a crystalline sample of disordered regions randomly 
distributed in space. The averaged picture for such a turbulent solid is exemplified by 
an exactly solvable lattice-gas model. The origin of heterophase turbulence is connected 
with stochastic instability of quasi-isolated systems. 



1 



1 Typical Stages of Motion 



The term turbulence is usually applied to characterizing chaotic motion of fluids [1]. There- 
fore, it could sound strange to mention turbulence in connection with solids. Nevertheless, 
a similar phenomenon can occur in the latter, being related to mesoscopic heterophase fluc- 
tuations [2]. 

We shall keep in mind everywhere below the word turbulence in its standard sense [3] as 
characterizing the motion possessing the following main stochastic properties: (i) sensitive 
dependence on initial conditions, due to exponential divergence of initially close trajectories; 
(ii) chaotic space structure at fixed moments of time; (iii) existence of a strange attractor. 
Note that Ruelle [4] used the title of turbulent crystal with regard to a hypothetical stable 
amorphous solid. However, real amorphous solids are always only metastable. 

Turbulence is a spatio-temporal phenomenon. In order to describe spatial and temporal 
scales typical of heterophase turbulence, let us recall what are, in general, characteristic 
scales for a system of interacting particles. 

The first evident characteristic length is an effective interaction radius r int or a scattering 
length a s ~ r int . If a characteristic particle velocity is v, then the interaction time is 



In the temporal interval < t < T int , called the dynamic stage, particles move, to a first 
glance, randomly. 

Another important length is a mean free path \ m f P defining the local- equilibrium time 



And the temporal diapason Ti nt < t < ti oc is the kinetic stage. At this stage, because 
of mutual interactions, interparticle correlations begin arising. For the density of particles 
p ~ a~ 3 , with a being a mean interparticle distance, the mean free path is A ~ (pa^) _1 , that 
is, A ~ a? /a 2 s . Therefore, 



In condensed matter, a ~ r int , because of which T ioc ~ r int , and the kinetic stage being 
practically absent. 

After the local equilibrium time Ti oc , the state of local equilibrium is being evolved, and for 
t > Tioci one has the hydrodynamic stage. But the latter consists of two qualitatively different 
parts. There exists a coherence radius r coh characterizing the effective size of mesoscopic 
heterophase fluctuations. The related coherence time is 



Heterophase fluctuations are mesoscopic because the coherence radius is between the mean 
interparticle distance a and the linear size of the whole system L, 



Tint 

Tint ~ 



V 




r coh 



7~coh 
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a < r coh < L . 
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The coherence time r co h ^> t\ oc defines an average lifetime of a mesoscopic fluctuation. The 
temporal interval Ti oc < t < r coh corresponds to the mesoscopic stage, when heterophase 
fluctuations look as frozen, making the system nonuniform. 

For times t > T co h, heterophase fluctuations become dynamic. They arise randomly in 
time, being chaotically located in space. Their appearance signifies the existence of local 
instabilities leading to divergence of initially close trajectories. At each fixed moment of 
time, the spatial structure of the system is chaotic. Dynamics of a system with heterophase 
fluctuations corresponds to the evolution on a strange attractor [2]. These characteristics 
are the main features of a turbulent motion [3] . This is why the occurrence of heterophase 
fluctuations can be called heterophase turbulence. This phenomenon can be noticed only on 
the macroscopic stage during times t > r coh . 

Heterophase fluctuations are rather ubiquitous in nature [2]. They can be clearly ob- 
served by various experimental techniques. For instance, in many magnetic alloys, they have 
been studied by neutron scattering [5]. Mesoscopic heterophase fluctuations are common 
in liquid crystals [6], high-temperature superconductors [7-10], and many other materials 
[2]. Anomalous atomic diffusion in solids [11] and anomalous transport in some dynamic 
systems [12] can be explained by mesoscopic fluctuations. Such an anomalous diffusion is an 
analog of the turbulent diffusion [3]. Highly unusual properties of colossal magnetoresistive 
manganites are due to the coexisting fluctuating domains of semiconducting nonmagnetic 
states and metallic ferromagnetic states [13]. A good example of such a colossal magne- 
toresistive manganite is Lax^Ca^MnOs, with x ~ 0.33, where the fraction of the material 
in each of two well-defined states strongly depends on temperature and magnetic field, the 
inhomogeneity being characterized by a mesoscopic scale of about 100 A. Dynamical Jahn- 
Teller distortions around the antiferromagnetic to paramagnetic phase transition in UO2 [14] 
and N p 02 [15] are also a kind of mesoscopic fluctuations. Such local fluctuations arise in 
many liquids, including water, as mesoscopic cooperative structures of different polymor- 
phism [2,16-20], making a system dynamically heterogeneous. Seismic precursors related 
to local polymorphic transformations [21] can be viewed as mesoscopic heterophase phe- 
nomena. This concerns as well the appearance of droplets of quark-gluon plasma in a bulk 
hadronic medium and of hadronic bubbles in a bulk quark-gluon plasma, occurring around 
the quark-hadron transition [22,23]. Let us stress that heterophase fluctuations should not be 
confused with formally introduced microscopic clusters in the so-called equilibrium cluster 
approximations [24,25], but these fluctuations represent nonequilibrium mesoscopic germs 
with properties qualitatively differing from the surrounding matter. 

2 Averaged Picture of Turbulence 

There exists a statistical approach to describing turbulent motion in terms of an averaged 
picture [1,3]. Similarly to this, it is possible to develop an averaged picture of heterophase 
turbulence by averaging over heterophase fluctuations [2]. 

The appearance of turbulence in fluids happens when the Reynolds number surpasses a 
critical value. Analogously, heterophase turbulence arises only under a specific combination 
of interaction parameters and external variables, such as temperature or pressure. The 
Reynolds number for fluids is the ratio vL/i>, in which v is a characteristic velocity of motion, 
L is a characteristic linear size, and v is kinematic viscosity. That is, the Reynolds number is 
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the ratio of factors favoring turbulence to those stabilizing laminar motion. A counterpart of 
the Reynolds number for heterophase turbulence could be the ratio of interaction parameters 
favouring a disordered state to those stabilizing an ordered state. Although the definition of 
such a typical ratio for heterophase turbulence is not always simple and, in addition, depends 
on thermodynamic variables. 

In crystals, heterophase turbulence is manifested by the emergence of chaotically dis- 
tributed regions of crystalline disorder [16,17]. Factors favoring the appearance of such 
regions of disorder could be a high concentration of defects and strong anharmonicity. Good 
candidates for this could be quantum crystals [26] and strongly anharmonic crystals [27,28]. 
The realistic description of such crystals involves rather complicated techniques. In order 
to avoid mathematical complications, we shall resort here to a simple lattice-gas model [29] 
serving as a cartoon of a crystal with many defects or of a solid with pores and cracks [30]. 
An advantage of employing this model is the possibility to suggest its rigorous analysis and 
to explicitly demonstrate that heterophase turbulence can correspond to a more stable state 
than a uniform homophase state. 

Consider a system of N particles in volume V. Particles are distributed on a crystalline 
lattice so that each lattice site can be occupied by not more than one particle. The number 
of lattice sites N s can be larger than that of particles, N s ^> N; thus vacant sites can be 
present. Assume that the distribution of vacancies is not uniform over the system but that 
vacancies fluctuationally group in some parts of the sample, as a result of which the system 
possesses essentially different density of particles in different spatial regions, the latter being 
distributed chaotically. Such a picture is typical of heterophase turbulence, when a more 
dense phase and more dilute phase randomly fluctuate in space and time. 

In an averaged description, the total number of particles can be presented as the sum 
N = Ni + N 2 of particles pertaining to the dense phase, Aq, and of those in the dilute phase, 
N 2 . The whole volume of the system consists of two parts, V = V\ + V 2 , those occupied 
by the dense and dilute phases, respectively. Each volume V u , with v — 1,2, does not form 
a connected part filled by the corresponding phase, but is composed of a large number of 
disconnected mesoscopic regions randomly located in space. 

Since the number of sites, N s , is larger than the total number of particles, N, the density 
of sites is higher than that of particles, 



By definition, the dense phase is characterized by a higher density of particles that the dilute 
phase, 



Ps> P , 




N 



(1) 



Pi>p 2 , Pu = TT (i/ = l,2). 



(2) 



V 



The geometric probabilities of the corresponding phases are 




(w ± + w 2 — 1 , < w u < 1) . 



(3) 



The average density of particles can be written as 



p = w x p x + w 2 p 2 . 



(4) 



4 



It is convenient to introduce the dimensionless densities 

P 



Pv 

n u = — 

Ps 



n 



Then Eq. (4) reads 



n = WiUi + U>2^2 • 



(5) 



(6) 



To an i-site of the lattice, pertaining to a //-phase, one ascribes the variable taking 
the value = 1 when the site is occupied by a particle and the value = if the site is 
empty. The phase densities can be presented through the averages 



Pv 



1 



N s 



< Civ > , 
i=l 



N. 



(7) 



s i=i 



where < . . . > implies a statistical averaging. 

The procedure of averaging over heterophase fluctuations [2], for the considered case, 
yields the effective heterophase Hamiltonian 



H = H 1 ®H 2 , 

being a direct sum of the phase-replica Hamiltonians 



N a 1 N s 

H v = w u Y^ (Ki - /i) e iv + - wl Y ®ije iu e ju , 



1=1 



(8) 



(9) 



in which Ki is kinetic energy, \x is chemical potential, and $y is interaction potential. The 
phase probabilities w v are defined by minimizing the grand potential 



n = -e lnTre"^ {{3Q = 1) , 
where is average temperature in energy units. With the notation 

w = w± , W2 = 1 — w , 

the minimization equations are 

d 2 Q 



dw 



> . 



The first of these equations leads to the equality 

- — > = 
dw 

The second of Eqs. (12) tells that the heterophase susceptibility 

Xe 



1 d 2 n 

N s dw 2 > 



(10) 
(11) 
(12) 

(13) 

(14) 



is to be positive. Differentiating the thermodynamic potential (10) gives 



Xe 



< 



d 2 H 

dw 2 



> -p< 



'dIP 
dw , 



> 



(15) 
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The inequality xe > is a condition of heterophase stability. 

In order to illustrate in an explicit way the thermodynamic properties of the consid- 
ered system, without plunging to mathematical complications, let us choose the long-range 
interaction potential <J>jj = $ij(N s ) depending on the number of lattice sites so that 



lim $^(Ag = 0, 



N s ^oo 



lim 

N s ^oo 



1 



< OO 



(16) 



This type of interaction potentials allows one to get, in the thermodynamic limit, asymptot- 
ically exact solutions [31], which is due to the asymptotically exact decoupling 



For the compactness of the following expressions, let us introduce the notation 



1 



k = — y Ki , $ = — Y $a , 



s i=l 



iV, 



(17) 



(18) 



and let us define the dimensionless temperature and an effective chemical potential, respec- 
tively, 

/I- K 

/i* = 



T = — 

Also, denote the dimensionless densities (5) as 

rii = a, n 2 = b (a > b) . 
Then for the grand potential (10), one obtains 

^ r a -(i-aw + b - { i-b ){ i- w y-l^ 



(19) 



(20) 



-Tin 



4cosh 



aw — fi*w 



cosh 



6(1 - w) 2 - ji^l - w) 



(21) 



2T J 2T 
The heterophase susceptibility (15) can be reduced to the form 

Xe = (nl + nj)$ = (a 2 + b 2 )$. (22) 

The condition of heterophase stability, xe > 0, is nothing but the simple inequality 

$ > . (23) 

Heterophase turbulence is possible for positive interaction strength $, which plays here the 
role of the Reynolds number. 

The condition of heterophase stability (23) is a necessary, but not sufficient, condition for 
the occurrence of heterophase turbulence. To show that the latter really develops, one has to 
check the existence of solutions for all thermodynamic characteristics. We have accomplished 
such an analysis by means of numerical calculations. Figure 1 demonstrates the temperature 
dependence of the partial densities (20) for different average densities n, and Fig. 2 presents 
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the temperature dependence of the probability of the dense phase. Physically meaningful 
solutions require that the probability w be in the interval < w < 1. Hence, heterophase 
turbulence can actually arise not for all densities and temperatures. The solutions for the 
partial densities (20), defined by the equations 



1 

2 



1 



tanh 



b(n — b) 



1 1 

- — - tanh 

2 2 



a(a — n 



have sense only when the dense-phase probability 

n — b 



w 



a — b 



lab — n(a + b) 
2(a-bfT 

2ab — n(a + b) 
2{a - bfT 



(a > b) 



(24) 



(25) 



satisfies the inequalities < w < 1. The latter, together with the condition a > b, yield the 
requirements 

b<n<a, 0<n< 0.323 . (26) 
Figures 3 and 4 show the effective chemical potential 



A** = 



\n — b) + b 2 {n — a) 
(a - b) 2 



and the dimensionless pressure 



P = - 



i dn 



p s $ dV 

as functions of temperature. The dimensionless specific heat 



and isothermal compressibility 



Cy = 



Kt = — 



T d 2 tt 



iv s $ or 2 




(27) 



(28) 



(29) 



(30) 



are shown in Figs. 5 and 6, respectively. Their positiveness demonstrates the thermal 
and mechanical stability of the heterophase system. Finally, Fig. 7 demonstrates that the 
dimensionless grand potential u w = Q/N s § for the heterophase matter is lower than that for 
the pure system, cui, when w — 1. This means that the state with heterophase turbulence is 
globally more stable than the state having no heterophase fluctuations. 

In this way, there exists a region of thermodynamic variables, e.g. temperature and 
density p, where the lattice system displaying heterophase turbulence is absolutely stable 
as a whole, despite that the emerging turbulence assumes local instability. That is, local 
instability renders the system more stable globally. 
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3 Origins of Heterophase Turbulence 



In general, there could be three causes producing heterophase turbulence, extrinsic, intrinsic, 
and one intermediate between the first two, which may be called stochastic [2]. 

The extrinsic cause presupposses the action of sufficiently strong external forces making 
the system nonequilibrium. These could be either specially designed time-dependent pertur- 
bations or uncontrolled random interactions with surrounding. This is similar to the way of 
preparing a turbulent fluid by strongly shaking the system [1], as a result of which the fluid 
becomes locally turbulent, i.e. nonuniform, but on the average it is uniform and isotropic, 
with zero average velocity. In the case of a solid, regions of disorder inside it could be formed 
by shock waves, irradiation by particle beams, electromagnetic fields etc. The possibility of 
producing heterophase turbulence by means of strong external forces looks rather evident, 
hence not as interesting. A more intriguing question is whether heterophase turbulence could 
emerge without strong external fields or not. 

The intrinsic mechanism for the formation of heterophase fluctuations is based on the 
Boltzmann hypothesis [32] that in the entirely equilibrium system strongly nonequilibrium 
local fluctuations can appear in a spontaneous way. This mechanism of the appearance of 
nonequilibrium heterophase fluctuations from an equilibrium state was also advocated by 
Frenkel [33]. However, there is a slight inconsistence in this way of thinking. Really, first, 
one presupposses the existence of equilibrium states. Then one assumes that an equilibrium 
state can spontaneously pass to a nonequilibrium one. But the latter implies that equilibrium 
states are not absolutely stable, that is, they, actually, do not exist as stable states. Following 
this way of reasoning, it would be more logical to postulate from the very beginning that 
an absolutely stable state, to which an isolated system tends with time, is, generally, not an 
equilibrium state but a quasiequilibrium heterophase state [2]. 

The stochastic source generating heterophase turbulence is intermediate between the 
above two in the following sense. It accepts that completely isolated systems do not exist 
in reality, but any given system is always subject to uncontrollable random perturbations 
from the external surrounding. At the same time, one can prepare an almost isolated, or 
quasi-isolated, system for which external stochastic influence, if not completely eliminated, 
but is essentially reduced. Nevertheless, even infinitesimally small random perturbations can 
provoke noticeable chaoticity in a quasi-isolated system, with chaotic features governed by 
the system internal properties [2,34]. 

The fact that no real system can be ideally isolated has been repeatedly emphasized 
[35-39]. Moreover, it has been stressed [40,41] that the concept of an isolated system is 
self-contradictory by its own - To realize the isolation, one has to employ technical devices 
acting on the system, and to check that the latter is kept isolated, one needs to resort to 
measuring instruments perturbing the system. 

In this section we show that even negligibly small stochastic perturbations can leed to 
the development of chaoticity in the system. The evolution equations are, generally, partial 
differential equations. Let a dynamical system be described by a set of y(t) = [yi(x,t)\ of 
functions yi(x, t) enumerated by % — 1, 2, . . . , d, with a continuous variable x G © denoting a 
set of spatial variables pertaining to a domain D, and t G R+ being time. The set y(t), called 
the dynamical state, can be treated as a vector-column labelled by the discrete index i and 
by the continuous multivariable x. The velocity field v(y,t) = [vi(x,y,t)] is also a column 
in the space of i and x. Any system of evolution equations can be presented in the normal 
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form 

— y(t) = v(y,t) . 

The velocity field v(y,t) may contain differential and integral operations. The appearance 
of chaos is related to the problem of stability of evolution equations. 

3.1 Autonomous systems 

Before passing to more complicated cases, let us recall just in a few lines how the problem of 
stability is treated for autonomous systems of nonlinear partial differential equations. This 
introductory consideration is useful as far as in the standard textbooks one usually studies 
only the stability of nonlinear ordinary differential equations. At the same time this example 
will show how to deal with the matrix notation [34] that will be constantly used in what 
follows. 

A dynamical system is autonomous if the velocity field does not depend explicitly on 
time, v(y,t) = v(y). Then the system evolution is given by a set of equations 

j f V(t) = v(y) , (31) 

presented in the matrix form. 

One looks for a stationary solutions y* defined by the zero velocity field, 

v(y*) = , y* = [y t (x)] . (32) 

For a small deviation 

5y*(t)=y(t)-y* (33) 
from the stationary solution y* , one gets a linearized equation 

j t 5y\t) = J*5y*(t) , (34) 

in which J* = [J*j(x, x')] is a contraction matrix with the elements 

J ^ x) = ^jW (35) 

Looking for the solutions of the eigenproblem 

J*fn = J* n fn , (36) 

where n is a multi-index, one finds the eigenvalues J* and the eigenfunctions f n = [f n i(x)] of 
the contraction matrix. Assume that the eigenvalues form an orthonormalized and complete 
basis, {f n }, so that 

fmfn 3mn j fnf n ^ ■ 

n 

Note that in the explicit presentation, the matrix eigenproblem (36) reads 

J J ij( X ' X ') fnjW) dx' = J*f ni (x) ■ 



9 



Expanding the deviations 

6y*(t)=J2 Sc n(t)fn 
n 

over the basis {f n }, one has the equation 

^ Sc n (t) = J*5c n (t) 



(37) 



(38) 



for the coefficients 5c n (t). This equation immediately gives 

Sc n (t) = 5c n (t ) exp {J*(t - t )} . (39) 
The stability of the stationary solutions is characterized by the Lyapunov exponents 



A, 



lim - In 



5c n (t) 



Sc n (t ) 



Re Jl . 



(40) 



The set {A„} is the Lyapunov spectrum. The functions /„ are termed normal modes. One 
tells that an n-mode is stable if A n < 0, neutral when \ n = 0, and unstable for A n > 0. The 
motion is chaotic if at least one mode is unstable. 



3.2 Stochastic systems 

Let us now turn to the problem of stability of stochastic differential equations. Assume that 
the system is subject to the action of stochastic fields £(£). We continue using the matrix 
notation, where the dynamical state y((,,t) = [yi(x, £,i)] and the velocity field t>(y,£,t) = 
[vi(x, y, £, t)] are the columns in % and x. The evolution of a stochastic dynamical system is 
given by the matrix equation 

j t y(U) = v(y,U)- (41) 
Let the evolution equation (41) be complimented by the initial condition 

2/(^,0) = 2/(0) , (42) 

where y(0) is a set of given functions yi(x). To solve a stochastic differential equation implies 
to find the solution 

y(0=«y(e,i)» (43) 

averaged over the stochastic fields £(£). We denote the stochastic averaging by the double 
angle brackets <C . . . 3>, reserving the single brackets < . . . > for statistical averaging, as in 
Sec. 2. 

Note that stochastic differential equations [42] can be defined either in the sense of Ito 
or in the sense of Stratonovich. We prefer to use the latter definition which involves simpler 
calculations, agreeing with the standard analysis, and which is better motivated physically 
[43]. Another possibility could be to present the stochastic fields as expansions, with 
random coefficients over smooth functions of time [44,45]. This method makes it possible to 
employ the standard analysis in differentiating and integrating over time. The final results 
of the expansion approach [44,45] coincide with those of the Stratonovich way. 
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The stability of the stochastic equation (41) cannot be analysed in the same manner as 
the stability of the autonomous equation (31), since the stationary solutions for the stochastic 
equation (41) are not defined. Another approach is necessary. Let us introduce the local 
stability factor 

Ht) . sup §<M (44) 
sy(o) \oy{0)\ 

characterizing the maximal deviation of the average trajectory at time t after an infinitesimal 
variation of initial conditions 5y(0) — > 0. The following classification for the local property 
of motion is evident: 

v(t) < 1 locally stable; 
u(t) = 1 locally neutral; 

u{t) > 1 locally unstable . (45) 

If the limit 

u(oo) = lim fit) 

exists, then one can tell that the motion is asymptotically stable for z/(oo) < 1, asymptotically 
neutral for z/(oo) = 1, and asymptotically unstable for z/(oo) > 1. The quantity 

A = lim - In u(t) 

is the largest Lyapunov exponent. 

The usage of the local stability factor (44) provides one with a much richer information 
on the motion than the largest Lyapunov exponent, since it is often important to study local 
in time stability but not only the asymptotic, as t — > oo, property of motion [34]. Many 
dynamical systems display a very complicated structure of their phase space, containing 
different singular zones, such as trapping islands, because of which the fine local properties 
of orbits play a leading role, while fairly rough characteristics such as the Lyapunov exponent 
are less important [46]. In addition, for stochastic differential equations, the divergence of 
trajectories is not necessarily exponential [47]. For example, if v(t) ~ t 13 , then the motion 
can be either asymptotically stable or unstable depending on whether (3 < or (3 > 0, while 
the Lyapunov exponent A = for both these cases characterizes the motion as neutral, thus, 
loosing information about the motion. 

To calculate the stability factor (44), we need to find 5y(t) =<^ 5y(£,t) ^$>. For this pur- 
pose, we introduce the stochastic multiplier matrix [34] as a matrix M(£, t) = [My(i, x', £, t)} 
with the elements 

^M'W$ (46) 

and the contraction matrix J(£, £) = [Jij(x,x',£,t)] with the elements 

J «( x > x lyJTJJ) ■ (47) 

Accomplishing the variational differentiation of the evolution equation (41), we get the equa- 
tion 

j t M(U) = J(t,t)M(t,t) (48) 
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for the multiplier matrix. When for the latter the eigenproblem 



M(£,t)<P n (t) = Vn(U)<Pn(t) ( 49 ) 

is defined, then the stability factor (44) becomes 

v(t) = sup| <A* n (£,0 > I • (50) 

n 

The eigenvalues fJ, n (£,t) of the stochastic multiplier matrix M(£,t) can be found from Eqs. 
(48) and (49). 



3.3 Quasi-isolated systems 

As is discussed in the beginning of this section, no real system can be ideally isolated, 
although the influence of random external perturbations could, probably, be made very small. 
Nevertheless, even arbitrarily small random perturbations could provoke system instability. 
To analyze this problem, we need to consider the stability of a given dynamical system 
not only with respect to the variation of initial conditions but also with respect to small 
stochastic perturbations. 

Let the amplitude of stochastic fields be small, which can be presented by writing an 
explicit small factor a in front of £(t). That is, instead of the evolution equation (36), we 
consider 

j£V(oit,t) = v(y,aZ,t) (a->0), (51) 

where the complex parameter a is infinitesimally small, so that a — > implies a,\ = Reo — > 
and o 2 = Ima — > 0. When a = 0, no stochastic fields act on the system, and the evolution 
equation 

j t y(0,t)=v(y,0,t) 

describes an isolated system. Switching on a small parameter models uncontrollable action 
of the surrounding. The parameter a is taken to be complex in order to simulate fluctuations 
of different physical quantities, for instance, of energy and attenuation, or of density and 
phase. 

The stability of a stochastic system is characterized by the stability factor (50), which 
for Eq. (51) takes the form 

v(e, t) = sup | < At n (a£, t) > | , (52) 

n 

where e = £(0:1,0:2) is a function of 01 = Reo and o 2 = Ima, such that e — > as a — > 0. 
Depending on how 01 — > and 0:2 — ► 0, the parameter e can behave either as e — > — or 
e — > +0, or it can even be that e = 0. Of all these possibilities, one has to choose the case 
corresponding to the maximal stability factor (52), as far as the latter, by definition, charac- 
terizes the maximal possible divergence of initially close trajectories. Therefore, everywhere 
in what follows, the limit e — > is defined as 

lim = sup lim . (53) 

£^0 e |e|->0 
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To analyze the asymptotic stability one has to look at the temporal limit t — > oo. The 
result of taking the double limit of the stability factor (52) may depend on the order in which 
the limits e ^ and t — > oo are taken. We shall say [34] that the system is stochastically 
stable if the order of these limits is not important, that is, the limits commute with each 
other, 

[lim ,lim]z/(e,t) = . (54) 
Conversely, the system is termed stochastically unstable, when these limits do not commute: 

[lim ,hm>M)^0. (55) 

If a quasi-isolated system is stochastically unstable then, since any real system is never 
completely isolated but can be only quasi-isolated, chaos in the system can be induced by 
infinitesimally small random perturbations. 



4 Stochastic Density Matrices 

Now, employing the technique of the previous section, let us investigate if the evolution of 
particle density could be stochastically unstable. If so, this would mean that the density 
as a function of space and time could exhibit nonequilibrium fluctuations. Each of such 
fluctuations, after appearing, lives a finite lifetime, and then disappears. But, as time tends 
to infinity, the density does not tend to an equilibrium, say uniform or strictly periodic in 
space, state. As far as the system is stochastically unstable, nonequilibrium fluctuations 
arise again and again, chaotically emerging in time as well as in space. The state of matter 
with such chaotic fluctuations would be representing heterophase turbulence. 
The density of particles 

p(r,t) = lim pi(r,r',t) , (56) 

r'^r 

as a function of the Cartesian vector r and time t, is a diagonal element of the single density 
matrix 

Pl ( r , r ',t)=<^( r ',t)^(r,t)> , (57) 

in which the brackets < . . . > imply statistical averaging and ip(r, t) is a field operator, being 
a vector-column in the space of internal degrees of freedom, such as spin, isospin, and so on. 
The evolution of the single density matrix (57) is coupled with that of the double density 
matrix 

p 2 ( ri , r 2 , ri, r' 2 , t) =< ^(r' 2 , t)^(r[, t)i/j(r 1: i)^(r 2 , t) > . (58) 

General properties of density matrices can be found in the book [48] (for some recent devel- 
opments see Refs. [49,50]). 

The evolution equation for the density matrix (57) can be presented [39] in the form 



1 Ql Pl( r l> r 2,£) 



— (V\-Vt) + U{v u t)-U{v 2 ,t) 



Pi(ri,r 2 ,t) + 



2m 

+ J [$(ri, r 3 ) - $(r 2 , r 3 )] p 2 (r 3 , n, r 3 , r 2 , t) dr 3 , (59) 

where we set the Planck constant H = 1, m is particle mass, and U(r,t) is an external 
potential, while $(r, r') = $(r',r) is an interaction potential. 
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When r 2 — > ri, then Eq. (59) reduces to the asymptotic form 

^ Pi(r 1? r 2 ,t) = ^- (V? - V*) Pl ( ri ,r 2 ,t) . 
Taking account of the equality 



(60) 



V\-V 2 2 = (Vi + V 2 ) (Vi - V 2 ) 
and introducing [39] the density-of-currrent matrix 



j(ri,r 2 ,*) 



2mt 



(Vi - V 2 )pi(ri,r 2 ,t) 



one can rewrite Eq. (60) as 

d_ 
di 



Pl (r 1 ,r 2 ,t) = -(Vi + V 2 )j(r 1 ,r 2 ,t) . 



The density of current is the diagonal element 

j(r,t) = lim j(r,r',t) , 

r'— >r 

which can also be represented as the statistical average 

j(r,i) =< j(r,t) > 

of the density-of-current operator 



(61) 



(62) 



(63) 



(64) 



j(r,0 = - 



2mr 



where, for brevity, we write ip = ip(r,t). In the limit r 2 = ri, Eq. (62) transforms to the 
standard continuity equation 



^ Pl (r,t)-V-j(r,t) = 0. 



(65) 



Hence, we may say that Eq. (60) asymptotically, as r 2 — > ri, equivalent to the continuity 
equaiton (65) defining the evolution of the particle density (56). 

To check the stochastic stability of the evolution equation (60), one has to incorporate 
there a small stochastic term. For the compactness of notation, let us denote x = {r!,r 2 }. 
Then an asymptotically small stochastic field is a£(x,t), with a — > 0. Also define 

y(x,a£,t) = pi(ri,r 2 ,a£,£) . (66) 

The stochastic equation, obtained from the evolution equation (60), reads 

d 



dt 



where the velocity field is 

v(x,y,a£,t) 



y(x,a£,t) = v(x,y,a£,t) , 



(V?-V!)+a£(z,f) 



2mi 



y{x,a£,t) . 



(67) 



(68) 
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The stochastic term initiates here fluctuations of energy and of attenuation leading to the 
related fluctuations of the amplitude and phase of the density matrix. 

Following the strategy of the previous section, we find the contraction matrix (47), 

J(x, x', < t) = (V? - Vl) + a£(x, f )] 5(x - x') . (69) 

Since a — > 0, the eigenvalues of the matrix (69) can be found by means of perturbation 
theory. In the zero approximation, the normalized eigenf unctions of the contraction matrix 
are 

<Pk(x) = y ex P W k i ' r i + k 2 • r 2 )} , (70) 

with the notation k = {k 1; k 2 }. The eigenvalues of the contraction matrix, in the first order, 
are 

J k (a£,t) = J cp* k (x)J(x,x',a^,t)cp k (x') dx dx' . (71) 
This, with the functions (70), gives 

JkH, t) = - ^- (k? - kfj + a£ k (t) , (72) 

where 

£ k (t) = / (pl(x)Z(x,t)(p k (x) dx . (73) 



Note that if £(x,t) = £(t) is a purely temporal perturbation, then Eqs. (70) and (72) give 
exact eigenfunctions and eigenvalues of the contraction matrix. 

The eigenfunctions (70) are stationary, which means [34] that in this approximation the 
multiplier matrix (46) possesses the same eigenfunctions. The eigenvalues of the multiplier 
matrix (46) are given by the eigenproblem 

J M(x, x\ a£, t)<pk(x') dx' = ti k (a£, t)<p k (x) . (74) 

Taking into account the equation (48) for the multiplier matrix, one has 

H k (a£,t) = exp j- ^- (k\ - kpj t + a j\ k {t') dt'} . (75) 

Assuming that £(x, t) is a real Gaussian stochastic variable, centered at zero, for the stability 
factor (52), having here the form 



v(e, t) = sup | < /ifc(«^, t) » | , 
k 



we obtain 



u(e, t) = sup exp |e J* < Ck(t')Ck(t") > dt' dt"} , (76) 

with e = \{pt\ — a|). Keeping in mind that 5^0 according to the limit (53), the factor 
(76) can be rewritten as 

i/(M)=exp{|e|r(i)i}, (77) 
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where 

r(t) = sup 1 

It is evident that 



JO 



(78) 



lim lim i/(e, t) = 1 . (79) 

However, under fixed \e\ > 0, it may happen that the limit of v(e,t) as t — > oo, does not 
exist, which could occur e.g. if T(t) oscillates. Or it may be that T(t)t — > oo, as i — > oo, 
which results in t) — > oo. The latter behaviour can be illustrated by assuming that 
£ (x, i) describes the temporal white noise with the correlation function 

<£(x,0£(^,O >= iS(t-if), 

where 7 > 0. Then T(t) = 7, and the factor (77) becomes 

u(e,t) = exp(\e\jt) . 
Another possibility could be to consider the uniform noise with the correlator 

<£OM)£(z',0 > = 7 2 ■ 
Then F(t) = 7 2 t, and the factor (77) reduces to 

u(e,t) =exp(|e| 7 ¥) . 

In any of such cases, 

lim lim v(e,t) = 00 , (80) 
which demonstrates stochastic instability of the system, in agreement with definition (55). 



5 Conclusion 

The considered simple cases show that there always exist such infinitessimaly small random 
perturbations which make the system stochastically unstable. In general, for stochastic 
instability, it is sufficient that 

lim F(t)t = 00 . 

There exists a wide class of noises, for instance, ranging between the white noise and the 
uniform noise, which could provide the validity of this limit. 

Instability of an evolution equation implies that there are no stationary solutions for this 
equation. A chaotic dynamical system can only possess a chaotic attractor. The above ar- 
guments show that the evolution of the particle density (56) can be stochastically unstable. 
Hence, the density cannot tend, with time, to a stationary limit describing a uniform density 
for liquids or gases and a periodic density for crystals. But the density p(r, t) will always 
depend on time, occasionally displaying random spatio-temporal fluctuations of nonequi- 
librium nature. To be noticeable as nonequilibrium, these fluctuations have to be at least 
mesoscopic. Such nonequilibrium mesoscopic fluctuations make the system heterophase. The 
state of matter with heterophase fluctuations represents heterophase turbulence. 
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In this sense, there should exist turbulent crystals, whose density displays chaotic spatio- 
temporal fluctuations corresponding to mesoscopic regions of disorder. A turbulent crystal, 
being observed during the times shorter or of the order of the characteristic lifetime of het- 
erophase fluctuations, r co h, would be seen as a matter, in which the spatial regions having 
approximately a crystalline structure with a type of lattice symmetry are randomly inter- 
mixed with disordered regions having no symmetry. The turbulent crystal, although being 
turbulent, is anyway a crystal because, being observed during the times much longer than 
the heterophase coherence time r co h, possesses on the average a periodic density, with a 
well-defined symmetry of crystalline lattice. The corresponding density, averaged over such 
long times t ^> r co h, can be written as 

p(r)= lim^ /%(r,0 dt' . 

t— >oo t JO 

Assuming heterophase quasiergodicity [2], the average density writes 

p(r) = wipi(r) + w 2 p 2 , 

where W\ and w 2 are statistical weights of the crystalline phase and of the disordered phase, 
respectively. The crystalline density pi(r) is periodic, while that of the disordered phase is 
uniform. For the densities averaged over space, 




we return to the equality (4). 

The possibility for the existence of turbulent crystals does not require that all crystals 
be necessarily turbulent. Whether they are such or not depends on the internal properties 
of the system as well as on external conditions, as is demonstrated in the model picture of 
section 2. But this possibility to become turbulent must be taken into account. Heterophase 
turbulence may be absent at some particular conditions, but may emerge when changing the 
latter. For example, it may not exist at low temperatures, but can arise near the temperature 
of melting [16,17], where the very existence of the heterophase turbulence is necessary for 
understanding the physics of melting and crystallization. 
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Figure Captions 

Fig. 1. The partial densities of the dense phase (solid curve) and of the dilute phase 
(dashed curve), for different average densities n, as functions of temperature. 



Fig. 


2. 


The probability of the dense phase versus temperature for different n. 


Fig. 


3. 


The effective chemical potential vs. temperature. 


Fig. 


4. 


The dimensionless pressure as a function of temperature. 


Fig. 


5. 


The dimensionless specific heat vs. temperature. 


Fig. 


6. 


The isothermal compressibility vs. temperature. 


Fig. 7. 

pure dense 


The dimensionless grand potential for the heterophase matter, u w , and for the 
system, cu 1 , vs. temperature. 
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